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ABSTRACT 

I show that there is a physical limit to the mass of a black hole, above which it can- 
not grow through luminous accretion of gas, and so cannot appear as a quasar or active 
galactic nucleus (AGN). The limit is Mmax & 5 x 10'° Me@ for typical parameters, but can 
reach Mmax & 2.7 x 10!! Me in extreme cases (e.g. maximal prograde spin). The largest 
black hole masses so far found are close to but below the limit. The Eddington luminosity 
~6.5 x 10* erg s~! corresponding to Mmax is remarkably close to the largest AGN bolometric 
luminosity so far observed. The mass and luminosity limits both rely on a reasonable but 
currently untestable hypothesis about AGN disc formation, so future observations of extreme 
supermassive black hole masses can therefore probe fundamental disc physics. Black holes can 
in principle grow their masses above M,,,x by non-luminous means such as mergers with other 
holes, but cannot become luminous accretors again. They might nevertheless be detectable in 
other ways, for example through gravitational lensing. I show further that black holes with 
masses ~M,y,x can probably grow above the values specified by the black-hole—host-galaxy 
scaling relations, in agreement with observation. 
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1 INTRODUCTION 


Astronomers generally agree that the centre of almost every 
galaxy contains a supermassive black hole (SMBH), with masses 
M ~ 10°-10!°Mo. The observed hole masses correlate tightly 
with large-scale properties of the host galaxy’s central bulge (see 
Kormendy & Ho 2013, for a recent review). This initially surprising 
connection arises because the gravitational potential energy released 
as a black hole grows offers the most efficient way of extracting en- 
ergy from ordinary matter (e.g. Frank, King & Raine 2002), and 
could potentially destabilize a host galaxy’s central bulge (King 
2003). The huge luminosities they produce make accreting black 
holes detectable as quasars and active galactic nuclei (AGN), and ex- 
ert mechanical feedback on their surroundings. This regulates their 
growth rates, and limits their masses to values specified by proper- 
ties of the host (for a recent review see King & Pounds 2015). As the 
host galaxies grow their masses, the black holes can grow further. 
An obvious question is whether there is a limit to this process, or 
whether a black hole can in principle reach any given mass, given 
a suitable host and enough time. The first attempt to answer this 
question was made by Natarajan & Treister (2009). They derived a 
limit by arguing that a self—gravitating accretion disc would blow 
itself away above a certain black hole mass of order 10°Mo, the 
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precise value depending on the properties of the host galaxy’s dark 
matter halo, in agreement with two lines of observational evidence. 

I study this question further here. In Section 2, I give a simple 
argument for a physical limit on any SMBH mass, above which it 
cannot form a disc and so grow by luminous accretion, together 
with some direct consequences. In Section 3, I consider the effects 
of changing some of the assumptions made in the simple argument 
of Section 2, and Section 4 is a discussion. 


2 SUPERMASSIVE BLACK HOLE GROWTH 


2.1 Disc accretion in AGN 


It has long been known that SMBH grow their masses mainly by 
luminous accretion of gas (Soltan 1982). Since gas within galaxies 
must have significant angular momentum, SMBH accretion must 
at any given time proceed largely through a disc (more probably, a 
series of discs of varying orientation). Any SMBH disc is likely to 
be self-gravitating outside some radius R,, ~ 0.01—0.1 pe (Collin- 
Souffrin & Dumont 1990; Shlosman, Begelman & Frank 1990; 
Huré et al. 1994). To see this, I note that the condition for stability 
against self-gravity for a gaseous disc can be expressed as 


6x2 
mGd 


(Toomre 1964), where c,, (2 and & are the local sound speed, 
orbital frequency and surface density, respectively, and G is the 
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gravitational constant. Vertical force balance in a disc requires 
c; = H&2, where H is the disc semithickness (e.g. Pringle 1981; 
Frank et al. 2002). Using this in (1) gives the stability requirement 
Q? M 
2nG  2nR3 
Here, 9 = & /2H is the mean density of the disc, and I have used the 
Keplerian relation Q = (GM/R°)'/”, with R the local disc radius, as 
appropriate for a thin disc around a black hole of mass M (Pringle 
1981; Frank et al. 2002) at the last step. The disc mass interior to 
radius R is My ~ 27tR* Hp, so (2) can be expressed in the well- 
known form 





p< (2) 


H 
Mas —M (3) 
R 


(Pringle 1981). 

For the parameters appropriate to thin discs around SMBH (I 
discuss this in more detail in Section 3 below) the disc aspect ratio 
obeys H/R ~ 10~? (cf. Collin-Souffrin & Dumont 1990; King et al. 
2008). The full disc equations then give 


Rog = 3 x 10a ng (L/Leaa) 7’ My”! cm, (4) 


where L is the accretion luminosity and Lgqq the Eddington lumi- 
nosity (cf. Collin-Souffrin & Dumont 1990; King & Pringle 2007). 
Here, a = 0.la91,7 =0.1901,m = M /Meaa are the standard vis- 
cosity parameter, the accretion efficiency and Eddington accretion 
ratio, respectively, and My = M/10* Moc. Gas cooling in the outer 
regions of these discs is fast enough that self-gravity is likely to 
lead to star formation rather than increased angular momentum 
transport (Shlosman & Begelman 1989; Collin & Zahn 1999). The 
very small aspect ratios H/R and masses M, of AGN discs mean that 
self-gravity appears first in modes with azimuthal wavenumber m 
~ R/H ~ 10°. These produce transient spiral waves which initially 
transport angular momentum (Anthony & Carlberg 1988; Lodato & 
Rice 2004, 2005). In a disc which is locally gravitationally unstable 
in this way, most of the gas initially at radii R > Rgo is likely either 
to form into stars, or to be expelled by those stars which do form, 
on a near-dynamical time-scale (cf. Shlosman & Begelman 1989). 

The estimate (4) is almost independent of parameters, and so must 
apply to almost every SMBH. Encouragingly, Ryo is only slightly 
smaller than the inner edge ~0.03 pc of the ring of young stars seen 
around the black hole in the centre of the Milky Way (Genzel et al. 
2003), strongly suggesting that the most recent accretion event on to 
the central SMBH formed an accretion disc within R,, and passed its 
angular momentum to the self-gravitating region further out which 
produced these stars. The small disc masses My < 10~°M expected 
from the self-gravity constraint offer an immediate explanation of 
the preferred time-scale ~10° yr of SMBH growth phases (King & 
Nixon 2015; Schawinsk1 et al. 2015). 

This reasoning implies that the outer radius of any SMBH ac- 
cretion disc cannot exceed Ry, which is effectively independent of 
the SMBH mass. But the inner disc radius must be at least as large 
as the ISCO Gnnermost stable circular orbit) around the SMBH, 
whose size scales directly with the SMBH mass M, as 


GM 3 
Risco = Ia = 3.7 310 Msg fs cm. (5) 


Here, f(a) is a dimensionless function of the SMBH spin parameter 
a, with f (a) = 5f5(a), so that fs ~ 1 corresponds to prograde accretion 
at moderate SMBH spin rates a ~ 0.6. If Risco Z Rsg, disc accretion 
is likely to be suppressed. Any disc material arriving at Risco feels 
only very weak outward angular momentum transport. If the SMBH 
accretes any of this matter at all, it must be self-gravitating and so 
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swallowed whole, without radiating as a disc. An SMBH might in 
principle grow its mass in this way, but we shall see below that it 
cannot subsequently reappear as a bright disc-accreting object, i.e. 
a quasar or AGN. 


2.2 SMBH mass limit 


Comparing (4) and (5), we see that the ISCO radius exceeds the 
self-gravity radius, making disc formation impossible, for SMBH 
masses larger than 


Mmax = 5 x 10" Moag, Nor (L/ Ley ——, (6) 


This is a physical upper limit to the mass of the SMBH in any quasar 
or AGN, since these systems have accretion discs. 

Fig. 1 shows the curve M = M,,,x(a), with a = 0.1 and L = Lgaa, 
while 7, f5(a) are specified parametrically as functions of a (cf. 
the relations 9 and 11 in King & Pringle 2006). The whole curve 
lies slightly above all the masses measured for accreting SMBH 
except for 0014 + 813 (Ghisellini et al. 2010; M ~ 4 x 10'°Mo) 
and H1821+643 (Walker et al. 2014; M ~ 3 x 10'° Mo). The first 
system is compatible with the limit provided that accretion is pro- 
grade and a = 0.2. H1821+643 is compatible with the limiting mass 
provided that a = —0.1 — that is, prograde accretion is possible for 
any spin parameter a > 0, but retrograde accretion on to this hole 
with |a| > 0.1 is ruled out. 

For spin rates a = | corresponding to maximal prograde spin 
wit the sense of accretion, the normalization in (6) becomes 2.7 x 
10'' Mo, which is the absolute maximum for an accreting SMBH. 
In practice values of M of this order are likely to be rather rare, as 
this requires disc accretion to be almost permanently prograde as 
the hole mass grows (the spin-down effect of retrograde accretion 
is greater than spin-up by prograde accretion, because of its larger 
ISCO and so its lever arm). This in turn probably requires the 
hole spin to be permanently correlated with a fixed direction of 
the potential controlling gas flow within the galaxy, and so would 
tend to produce a spin axis and hence AGN jet direction which is 
similarly aligned with the galaxy. Observations do not support this 
predicted correlation (Nagar & Wilson 1999; Kinney et al. 2000; 
Sajina et al. 2007). 

If accretion is not controlled by a large-scale potential in this 
way, it presumably has to involve multiple small-scale events, es- 
sentially random in time and orientation. This ‘chaotic’ type of 
accretion (King & Pringle 2006, 2007; King & Nixon 2015) leads 
statistically to spin-down, again because retrograde events have 
larger lever arms than prograde, and occur almost as often. King 
et al. (2008) show that this type of feeding predicts an attractor 
la} — ar 0.20M,,** for large SMBH masses (shown in Fig. 1). 
Other interactions with the SMBH also tend to reduce |a|. In partic- 
ular, mergers with other black holes statistically decrease the spin as 
|a| «x M~?* (Hughes & Blandford 2003). All these considerations 
suggest that SMBH usually cross the critical M = Mj. curve (6) 
at modest values of |a|, so that Minax & 5 x 10'° Mg in all but rare 
cases. 


2.3 AGN luminosity limit 


Assuming that the maximum observable luminosity of an AGN 
obeys the Eddington limit, the mass limit (6) implies a luminosity 
limit 


bing = 6:0 x 108 034? ng a ergs |, (7) 
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Figure 1. The mass limit Mmax for accreting supermassive black holes, compared with the largest observed masses. The curve shows Mmax as a function 
of black hole spin parameter a, where values a < 0 denote retrograde accretion. Accreting SMBH must lie below the curve. In order of decreasing mass, 
the systems shown are 0014 + 813, with M = 4 x 10/9 Mg (Ghisellini et al. 2010), the central quasar of the H1821 + 643 cluster (M = 3 x 10/0 Mo, 
Walker et al. 2014), NGC 4889 (McConnell et al. 2011), and the central galaxy of the Phoenix cluster (McDonald et al. 2012), both of these systems having 
M =2~ 10! Mo. Only the first two systems place any restrictions on the black hole spin. For 0014 + 813, accretion must be prograde, with a = 0.2. For 
H1821 + 643, if accretion is prograde |a| can have any value, but if accretion is retrograde |a| must be less than about 0.1. All the other systems are compatible 
with any spin value. The dotted curves show the statistical effect |a]| « M~*+ (Hughes & Blandford 2003) of black hole mergers on mass growth. The attractor 
la] —> a ~ 0.20M co for chaotic gas accretion is shown (dotted track) with the grey surround indicating the typical spread in |a| (King, Pringle & Hofmann 


2008). 


This prediction is in remarkably good agreement with the highest 
QSO luminosity found in the recent WISE survey of hot, dust- 
obscured galaxies (‘Hot DOGs’; Assef et al. 2015, fig. 4), which 
appears to be the largest AGN luminosity so far found. (The SEDs of 
Hot DOGs are generally dominated by a luminous obscured AGN.) 


3 RAISING THE LIMIT 


The mass and luminosity limits (6, 7) agree well with current ob- 
servations. Their derivations given above are simple enough that it 
is worth asking about the effect of varying some of the assumptions 
made there. 

First, we should note that (6, 7) are the observable limits for lumi- 
nous accretion at the luminosity L, not absolute limits on the black 
hole mass (as I have remarked already, non-luminous mass growth 
beyond Minax 1s perfectly possible). Sub-luminous mass growth (i.e. 
at lower accretion rates, with L/Leaq < 1) beyond M,,,x 1s also pos- 
sible. Such systems would by definition be fainter and so harder 
to detect. A flare taking the luminosity briefly up to Lega does not 
help, as the SMBH mass would then exceed the M,,,, value corre- 
sponding to this brighter state and so not be detectable (physically, 
the higher accretion rate moves the self-gravity radius inside the 
ISCO, preventing luminous accretion). In a very similar way, other 
sub-luminous accretion, e.g. via an ADAF, does not get around the 
observable limit. 


The most radical way around the limits (6, 7) is fairly obvious on 
looking at equation (3). Luminous disc accretion with a significantly 
larger scaleheight H than the very small values (WH ~ 10~°R) for 
standard thin-disc AGN accretion could increase M,,,x dramatically. 
The most likely way this could potentially occur is if radiation 
pressure would dominate gas pressure in the outer parts of the disc, 
i.e. if much of the outer disc would lie in the region ‘a’ of Shakura 
and Sunyaev, rather than region ‘b’ (gas pressure dominant) as 
assumed above. The discussion by Kawaguchi, Pierens & Huré 
(2004) shows that this might happen near the self-gravity radius if 
the prevailing accretion rate satisfies 


M > Mya = 4a? Me yr! (8) 


which is sub-Eddington for SMBH masses M > 4 x 10%a}'} Mo. 
But is well known that disc region ‘a’ is strongly unstable on 
a thermal time-scale, both in the context of the a-prescription 
(Lightman & Eardley 1974) and in shearing-box simulations (Jiang, 
Stone & Davis 2013). Theory and simulations are not currently able 
to work out the consequences of this instability, so we should ask 
what observational constraints exist. If region ‘a’ is fed material 
from a disc region ‘b’ outside it, it evidently finds a way to supply 
matter to the black hole. We observe Eddington-limited systems at 
a wide range of black hole masses, supplied by a stable and long- 
lasting reservoir such as a companion star in a stellar-mass binary, 
sO an inner region ‘a’ is perfectly compatible with AGN feeding. 
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But it is much less obvious that an SMBH disc can form with region 
‘a’ conditions at its outer, self-gravitating radius, and stably feed 
the hole. An AGN feeding event probably involves a ballistic flyby 
of a mass of gas, dust and possibly stars, which becomes bound 
to the SMBH because of internal dissipation as tides act on it (cf. 
King & Pringle 2006; King & Nixon 2015). With radiation pressure 
already dominant at the self-gravity radius this dissipation seems 
likely to drive mass off rather than produce efficient AGN feeding. 
If accordingly we assume that the outer parts of discs feeding AGN 
cannot be in region ‘a’ conditions we are again left with the limits 
(6, 7). 


4 DISCUSSION 


There are several points to note about the limits Max, Lax. 

1. Once M > Myx the SMBH can still go on growing its mass, as 
I have remarked above, provided this does not involve luminous disc 
accretion. Indeed gas accretion by swallowing stars would produce 
very little radiation, removing the Eddington limit as a barrier to 
growth. But the arguments above suggest that this mass growth is 
unlikely to allow the black hole to reappear as an accreting quasar. 
An increase in |a| is unlikely, and disc accretion is in any case 
impossible for M > 2.7 x 10'' Mo@ for any value of a. One might 
nevertheless detect SMBH with masses above M,,,x 1n other ways, 
perhaps through gravitational lensing for example. 

2. At masses close to but below M,,,,, luminous disc accretion in 
a field galaxy is likely to approach the Eddington luminosity only 
rarely, since even the dynamical infall rate Lo /G (with f, the gas 
fraction) is below the Eddington value except in galaxy bulges with 
very high velocity dispersions o > 400kms~!. Even for bright- 
est cluster galaxies in the centres of clusters, accretion of cluster 
gas may be vigorous, but strongly super-Eddington rates appear un- 
likely. The sub-Eddington accretion likely to prevail in such systems 
does not trigger the strong feedback which probably underlies the 
M-o relation, so SMBH close to M,,,, can evolve above M—o, and 
need not make their hosts red and dead. The observational data (e.g. 
McConnell et al. 2011) suggest that this is indeed what happens. 

3. In line with this, M,,,.x lies well above the M-—o relation, since 
host bulge velocity dispersions do not reach the required values o 
> 700kms7!. 

4. The limits Myx, Emax might be breached if AGN discs could 
form with radiation pressure dominant at the self-gravity radius. So 
the survival of this limit in the face of further observations may 
have something to tell us about AGN disc formation, and indeed 
the non-linear development of the radiation pressure instability in 
these discs. 
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